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Current-induced dynamics in spin valves including composite free layer with antiferromagnetic 
interlayer exchange coupling is studied theoretically within the diffusive transport regime. We show 
that current-induced dynamics of a synthetic antiferromagnet is significantly different from dynam- 
ics of a synthetic ferrimagnet. From macrospin simulations we obtain conditions for switching the 
composite free layer, as well as for appearance of various self-sustained dynamical modes. Numer- 
ical simulations are compared with simple analytical models of critical current based on linearized 
Landau-Lifshitz-Gilbert equation. 
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I. INTRODUCTION 

After the effect of spin transfer torque (STT) in 
thin magnetic films had been predicted^ 2 - and then ex- 
perimentally proveniM, it was generally believed that 
current-controlled spin valve devices would replace soon 
the memory cells operated by external magnetic field. 
Such a technological progress, if realized, would certainly 
offer higher data storage density and faster manipula- 
tion with the information stored on a hard drive memory. 
However, it became clear soon that some important is- 
sues must be solved before devices based on spin torque 
could be used in practice. The most important is the 
reduction of current density needed for magnetic excita- 
tion (switching) in thin films, as well as enhancement of 
switching efficiency and thermal stability. Some progress 
has been made by using more complex spin valve struc- 
tures and/or various subtle switching schemes based on 
optimized current and field pulses^— . 

A significant enhancement of thermal stability can be 
achieved by replacing a simple free layer (single homoge- 
neous layer) with a system of two magnetic films sepa- 
rated by a thin nonmagnetic spacer, known as composite 
free layer (CFL). The spacer layer is usually thin so there 
is a strong RKKY exchange coupling between magnetic 
layer o 9 i 10 . In practice, antiferromagnetic configuration 
is preferred as it reduces the overall magnetic moment 
of the CFL structure and makes the system less vulner- 
able to external magnetic fields and thermal agitation. 
When the antiferromagnetically coupled magnetic layers 
are identical, we call the structure synthetic antiferro- 
magnet (SyAF). If they are different, then the CFL has 
uncompensated magnetic moment and such a system will 
be referred to as synthetic ferrimagnet (SyF). 

Current and/or field induced dynamics of CFLs is cur- 
rently a subject of both experimental and theoretical in- 
vestigations^—. Switching scheme of SyAF by mag- 
netic field pulses has been proposed in Ref. PJ|, and 
then the possibility of current-induced switching of SyAF 
was demonstrated expcrimcntally-i2. In turn, the possi- 



bility of critical current reduction has been shown for a 
CFL with ferromagnetically coupled magnetic layers^. 
However, the reduction of critical current in the case 
of antiferromagnetically coupled CFLs still remains an 
open problem. In a recent numerical study on switching 
a SyAF free layer— it has been shown that the corre- 
sponding critical current in most cases is higher than the 
current required for switching of a simple free layer, and 
only in a narrow range of relevant parameters (exchange 
coupling, layer thickness, etc.) the critical current is re- 
duced. Hence, proper understanding of current-induced 
dynamics of CFLs is highly required. We also note, that 
CFL can be used as a polarizer, too. Indeed, it has been 
shown recently^ that SyAF used as a reference layer 
(with magnetic moment fixed fixed to adjacent antifer- 
romangentic layer due to exchange anisotropy) might be 
excited due to dynamical coupling^ with a simple sens- 
ing laye r 19 : 20 . 

The main objective of this paper is to study current- 
induced dynamics of a CFL with antiferomagnetic 
RKKY coupling in metallic spin valve pillars. We con- 
sider a system AF/F /Ni/Fi/N 2 /F 2 , shown in Fig. [TJ 
where AF is an antiferromagnetic layer (used to bias 
magnetization of the reference magnetic layer Fo), Fi 
and F 2 are two magnetic layers, while Ni and N2 are 
non- magnetic spacers. The part Fi/N 2 /F 2 constitutes 
the CFL structure with antiferromagnetic interlayer ex- 
change coupling. We examine current-induced dynamics 
of both SyAF and SyF free layers. These two structures 
differ only in the thickness of Fi layer, while RKKY cou- 
pling and other pillar parameters remain the same. 

We assume that spin-dependent electron transport is 
diffusive in nature, and employ the model described in 
Refs. [2l| and [22]. An important advantage of this model 
is the fact that it enables calculating spin current com- 
ponents and spin accumulation consistently in all mag- 
netic and nonmagnetic layers, as well as current-induced 
torques exerted on all magnetic components. The torques 
acting at the internal interfaces of CFL introduce ad- 
ditional dynamical coupling between the corresponding 
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FIG. 1. Scheme of the spin valve pillar structure with a com- 
posite free layer. 



magnetic layers. Consequently, the magnetic dynamics 
of CFL has been modelled by two coupled macrospins 
and described in terms of Landau-Lifshitz-Gilbcrt (LLG) 
equation. In addition, we derive some analytical expres- 
sions for critical currents from the stability conditions 
of linearized LLG equation in the static point o 23 ' 24 , and 
discuss results in the context of numerical simulations. 

The paper is organized as follows. In section 2 we de- 
scribe the assumed models for spin dynamics and STT 
calculations. In section 3 we analyze STT acting on 
CFL and present results from numerical simulations on 
current-induced switching and magnetic dynamics. Some 
additional information on STT calculation can be found 
in the Appendix. Critical currents are derived and dis- 
cussed in section 4. Finally, summary and general con- 
clusions are in section 5. 



II. MAGNETIZATION DYNAMICS 

In the macrospin approximation, magnetization dy- 
namics of the CFL is described by two coupled Landau- 
Lifshitz-Gilbert (LLG) equations, 
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for i = 1,2, where Si stands for a unit vector along the 
net spin moment of the i-th layer, whereas H c fa and Tj 
are the effective field and current-induced torque, respec- 
tively, both acting on Si. The damping parameter a and 
the saturation magnetization M s are assumed the same 
for both magnetic components of the CFL. Furthermore, 
7 g is the gyromagnetic ratio, (io is the vacuum perme- 
ability, and di stands for thickness of the F, layer. 
The effective magnetic field for the Fj layer is 
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where i, j = 1,2 and i =/= j. In the latter equation, i? app 
is the external magnetic field applied along the easy axis 
in the layers' plane (and oriented opposite to the axis 
z ), -ffani is the uniaxial anisotropy field (the same for 



both magnetic layers), and Hd cmt = (Nj • Si)M s is the 
self-demagnetization field of the Fj layer with Nj being 
the corresponding demagnetization tensor. Similarly, 
Hint; = (N 0i • S )M s0 + (Nji ■ Sj)M s describes the mag- 
netostatic influence of the layers Fo and Fj on the layer 
Fj, respectively. Here, M s q is the saturation magnetiza- 
tion of the layer Fo, which might be generally different 
from M s . Components of the tensors Nj, Noi and Njj 
used in our simulations have been determined by the 
numerical method introduced by Newell et al^Z. This 
method was originally developed for magnetic systems 
with non-uniform magnetization. To implement it into 
a macrospin model we considered discretized magnetic 
layers with uniform magnetizations, calculated tensors 
in each cell of the layer, and then averaged them along 
the whole layer. Since these tensors are diagonal, the de- 
magnetization and magnetostatic fields can be expressed 
a s Hdomi = (H^Six, HfySiy, Hf z Si Z ), and Hi n tj = 
(H x l Sox,H y l So y ,H z z Soz) + (H^Sjx, H^Sjy, H z l Sj Z ), 
with Six, Si v , and Si Z denoting the components of the 
vector Si (i = 0, 1, 2) in the coordinate system shown in 
Fig.l. Finally, -ffrcKKY i stands for the RKKY exchange 
field acting on Si, which is related to the RKKY 
coupling constant as H KKKYi = - Jrkky/CmoM^)^- 

To include thermal effects we add to the effec- 
tive field © a stochastic thermal field H t hi = 
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For both spins its compo- 



nents obey the rules for Gaussian random processes 
(H thi( (t)) = and (H thiC (t)H thje {t')) = 2DS tJ 5 a S(t - 
if), where i,j = 1,2 and f, £ = x,y,z. Here, D is the 
noise amplitude, which is related to the effective temper- 
ature, T ff, as 



D = 



(3) 



where fee is the Boltzmann constant, and Vj is the volume 
of Fj layer. 

In general, the current-induced torques acting on Si 
and S2 can be expressed as a sum of their in-plane and 
out-of-plane components t\ = t"i|| + tij_ and T2 =1^11 + 
T2j_, respectively. In a CFL structure, the layer Fi is 
influenced by STT induced by the polarizer Fo, as well 
as by STT due to the layer F 2 . In turn, the layer F 2 is 
influenced by the torques from the layer Fi. Hence we 
can write 
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where / is the charge current density, which is positive 
when electrons flow from the layer F2 towards Fo (see 

Fig- H]), while the parameters a?' and = 1,2) 

are independent of current /, but generally depend on 
magnetic configuration. 
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We write the current density in the spin space as j = 
jol + j • <r, where jo is the particle current density (7 = 
ejo), j is the spin current density (in the units of ft/2), 
a is the vector of Pauli matrices, and 1 is a 2 x 2 unit 
matrix. In frame of the diffusive transport model 2 ^, the 
parameters aj and fej are given by the formulas 



,(°) 
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where the angles # i an d #12 are given by cos # i = So ■ Si 
and cos 0i2 = Si • S2. Here, j' ly and are transver- 
sal (to Si) components of spin current in the layer Ni 
(taken at the N1/F1 interface) written in the local co- 
ordinate system of Si. Thus, j[ y is parallel to the vec- 
tor Si x (Si x So) (lying in the plane defined by So 
and Si) and j[ x is aligned with Si x So (perpendicular 
to the plane defined by So and Si). Note, that the z- 
component of spin current is aligned along Si and does 
not contribute to STT—. Similarly, we define the torque 
amplitudes acting inside the CFL. Here, 2% an d j'ix are 
the components of spin current in the layer N2 (taken 
at the N2/F1 interface), transversal to Si and lying re- 
spectively in the plane and perpendicularly to the plane 
defined by Si and S2. Analogically, j%y and j 2 " are spin 

current components in N2 transversal to S2 . The relevant 
spin current transformations are given in Appendix [X] 

Note, that the spin current components depend linearly 

(1) (i) 

on I, so the parameters a\ and b/ are independent of 
the current density I. 



III. NUMERICAL SIMULATIONS 

In this section we present results on our numer- 
ical simulations of current-induced dynamics for two 
metallic pillar structures including CFL with anti- 
ferromagnetic interlayer exchange coupling. As de- 
scribed in the introduction, the considered pillars 
have the general structure AF/Fo/Ni/Fi/N 2 /F 2 (see 
Fig.(l)). More specifically, we consider spin valves 
Cu - IrMn(10)/Py(8)/Cu(8)/Co(di)/Ru(l)/Co(d 2 ) - Cu, 
where the numbers in brackets stand for the layer thick- 
nesses in nanometers. The layer Py(8) is the Permalloy 
polarizing layer with its magnetization fixed due to ex- 
change coupling to IrMn. In turn, Co(di)/Ru(l)/Co(cfe) 
is the CFL (F1/N2/F2 structure) with antiferromagnetic 
RKKY exchange coupling via the thin ruthenium layer. 
The coupling constant between Co layers has been as- 
sumed as Jrkky — — 0.6mJ/m 2 , which is close to ex- 
perimentally observed value o 19 i 20 . Here, we shall analyze 
two different geometries of CFL. The first one is a SyAF 



structure with d\ = di — 2 nm, while the second one is a 
synthetic ferrimagnet (SyF) with d\ = 2d 2 = 4nm. 

Simulations have been based on numerical integration 
of the two coupled LLG equations (TT]) with simultaneous 
calculations of STT, see Eq.@. We have assumed typ- 
ical values of the relevant parameters, i.e., the damping 
parameter has been set to a = 0.01, while the uniax- 
ial anisotropy field i7 an i = 45kAm~ 1 in both magnetic 
layers of the CFL. In turn, saturation magnetization of 
cobalt has been assumed as M s (Co) = 1.42 x 10 6 Am -1 , 
and for permalloy M 8 (Py) = 6.92 x 10 s Am" 1 . The 
demagnetization field and magnetostatic interaction of 
magnetic layers have been calculated for layers of ellip- 
tical cross-section, with the major and minor axes equal 
to 130 nm and 60 nm, respectively. 

For both structures under consideration we have an- 
alyzed the current-induced dynamics as a function of 
current density and external magnetic field. The results 
have been presented in the form of diagrams displaying 
time-averaged values of the pillar resistance. Numerical 
integration of Eq.([T]) has been performed using corrector- 
predictor Heun scheme, and the results have been veri- 
fied for integration steps in the range from 10 -4 ns up 
to 10~ 6 ns. The STT components acting on CFL spins 
have been calculated at each integration step from the 
spin currents, which have been numerically calculated 
from the appropriate boundary conditions^. Similarly, 
resistance of the studied pillars has been calculated from 
spin accumulation in frame of the model used also for the 
STT description (for details see also Ref. H3) ■ 



A. Spin transfer torque 

Let us analyze first the angular dependence of STT 
components in the structures under consideration. Al- 
though the thicknesses of magnetic layers in the studied 
SyAF and SyF structures are different, the angular de- 
pendence of STT components as well as their amplitudes 
arc very similar. Thus, the analysis of STT in SyAF 
applies also qualitatively to the studied SyF free layer. 

First, we analyze STT components in the case when 
SyAF is rotated as a rigid structure, i.e. the antiparal- 
lel configuration of Si and S2 is maintained. To have a 
nonzero torque between Fi and F 2 layers, S2 has been 
tilted away from the antiparallel configuration by an an- 
gle of 1°. Figure [5] shows all three cartesian components 
(see Fig.l for definition of the coordinate system) of STT 
acting at N/F interfaces as a function of the angle 6 be- 
tween e z and Si. While the y and z-components are in 
the plane of the layers (the spins of CFL are rotated in 
the layer plane), the component x is normal to the layer 
plane. However, t x remains negligible at all interfaces of 
the CFL. The STT acting at N1/F1 reveals a standard 
(non-wavy^) angular dependence, and vanishes when Si 
is collinear with So . Its amplitude is comparable to STT 
in standard spin valves with a simple free layer. The STT 
at F1/N2 and N2/F2 interfaces also depends on the an- 
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FIG. 2. (color online) Angular dependence of the cartesian 
components of STT, in the units of hl/\e\, acting at N/F in- 
terfaces, when magnetization of the SyAF structure is rotated 
rigidly with both Si and S2 remaining in the corresponding 
layer planes. Here, 9 is an angle between Si and e z . S2 is 
tilted away from the antiparallel configuration with Si by an 
angle of 1°. 
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FIG. 3. (color online) Angular dependence of the cartesian 
STT components in the units of ft//|e| acting on F1/N2 and 
N2/F1 interfaces when So = Si — e z and S2 is rotated from 
e 2 in the layer's plane by an angle 6'. 

gle 8. However, they are about two order of magnitude 
smaller, which is a consequence of a small angle (1°) as- 
sumed between Si and S 2 . 

As will be shown in the following, CFL is usually not 
switched as a rigid structure, but generally forms a con- 
figuration which deviates from the antiparallel one. Fig- 
ure [3] shows how the STT components at the F1/N2 and 
N2/F2 interfaces vary when S2 is rotated from e z by 
an angle (?', while Si remains fixed and is parallel to 
So = e z . In such a case, the torque acting at N1/F1 
interface remains zero, as Si stays collinear to So- As 
before, the out-of-planc components are also negligible 
in comparison to the in-plane ones. The in-plane compo- 
nents of STT reveal standard angular dependence at both 
interfaces. The amplitude of STT at the internal inter- 
faces of CFL is comparable to that acting at the N1/F1 
interface in the case of noncollinear configuration of So 
and Si, when S2 is fixed in the direction antiparallel to 

When Si is noncollinear to So, the spin accumulation 



in Ni layer increases and consequently the amplitude of 
STT at Fi/N 2 and N 2 /F 2 decreases. In turn, when Si is 
antiparallel to So, the STT inside the CFL structure is 
reduced by more than a factor of 2. Nonetheless, the STT 
acting at the internal interfaces of the studied CFL layers 
might have a significant effect on their current-induced 
dynamics and switching process, provided the magnetic 
configuration of CFL might deviate remarkably from its 
initial antiparallel configuration. 



B. Synthetic antiferromagnet 

First, we examine dynamics of the SyAF free layer. 
From symmetry we have -Hrkky 1 = H RKKY 2 = ^rkky, 
and we have set Frkky = 2 kOe, which corresponds to 
Jrkky ~ — 0.6mJ/m . We have performed a number 
of independent numerical simulations modelling SyAF 
dynamics induced by constant current and constant in- 
plane external magnetic field. The latter is assumed to 
be smaller than the critical field for transition to spin- 
flop phase of SyAF. Accordingly, each simulation started 
from an initial state close to Si = — S2 = — e z . To have 
a non-zero initial STT for Si, both spins of the SyAF 
have been tilted by 1° in the layer plane so that they 
remained collinear. 

From the results of numerical simulations wc have con- 
structed a map of time-averaged resistance, shown in 
Fig. HJa). The resistance has been averaged in the time 
interval of 30 ns following initial 50 ns equilibration time 
of the dynamics. The diagram shows only that part of the 
resistance, which depends on magnetic configuration, and 
hence varies with CFL dynamics 2 ^. The constant part of 
resistance, due to bulk and interfacial resistances of the 
studied structure, has been calculated to be as large as 
R sp = 19.74 film 2 . For the assumed initial configuration, 
magnetic dynamics has been observed only for negative 
current density. When the current is small, no dynamics 
is observed since the spin motion is damped into the clos- 
est collinear state (Si = — S2 = — e 2 , marked as of 
high resistance. After exceeding a certain threshold value 
of current density, there is a drop in the averaged resis- 
tance, which indicates current induced dynamics of the 
SyAF free layer. Figurcs[4jb) and (f) show that this drop 
is associated with switching of the whole SyAF structure 
into an oposite state (Si = — S2 = e 2 , marked as t-J-) 

From Fig. 4(a) follows that the threshold current for dy- 
namics onset markedly depends on the applied field and 
reaches maximum at a certain value of H app7 H app = Ho- 
Furthermore, it appears that mechanisms of the switch- 
ing process for H app < Ho and H app > Hq are qualita- 
tively different. To distinguish these two mechanisms, we 
present in Figs. [4jb - i) basic characteristics of switch- 
ing, calculated for / = —1.0 x 10 8 Acm~ 2 and for fields 
H app = -400 Oe, which is below H [Figs, gfb - e)], 
and H app = 400, Oe, which lies above H [Figs. [4jf 
- i)]. Figs. BJb) and (f) present time evolution of 
the z-components of both spins. To better understand 
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FIG. 4. (color online) (a) Averaged resistance of Cu- 
IrMn(10)/Py(8)/Cu(8)/Co(2)/Cu(l)/Co(2)-Cu spin valve 
pillar with a SAF free layer as a function of current density 
and applied magnetic field. Examples of switching processes 
at I = -1.0 x 10 8 Acm~ 2 and # app = -400 Oe (b - e) and 
#a PP = 400 Oe (f - i). Figures (b) and (f) show dynamics 
of z-components of both spin moments, (c) and (g) present 
the overal magnetization of the free layer, (d) and (h) show 
the corresponding variation of pillar resistance, (e) and (i) 
show the spin trajectories of Si (red solid line) and S2 (black 
dashed line) in the time interval from to 10 ns, where switch- 
ing takes place. 

the SyAF dynamics, in Figs. SJc) and (g) we plotted 
the amplitude of overal SyAF magnetization, defined as 
in = \S\_ + S2 1 . This parameter vanishes for antiparallel 
alignment of both spins of CFL, but becomes nonzero 
when the configuration deviates from the antiparallel 
one. Magnetization of SyAF is also a measure of the 
CFL coupling to external magnetic field. Futhcrmorc, 
Figs.HJd) and (h) show the corresponding time variation 
of the resistance, 77, which might be directly extracted 
from experimental measurements as well. In addition, in 
Figs.QJe) and (i) we show the trajectories of Si and S2 in 
the real space taken from the time interval from t = to 
10 ns. In addition, from Figs.@|a) it has been found that 
the point where the threshold current reaches its maxi- 



mum is located at Hq ~ H ® 2 , which indicates its relation 
to magnetostatic interaction of F2 and fixed polarizer. 
This also has been confirmed by analogical simulations 
disregarding the magnetostatic coupling between mag- 
netic layers, which resulted in similar diagram, but with 
Hq = (not shown). This fact significantly facilitates 
understanding the mechanism of SyAF switching. 

The initial configuration assumed above was — Si = 
S2 — So with So = e z (4_t). When the magnitude of 
current density is large enough and 7 < 0, orientation of 
Si becomes unstable and Si starts to precess with small 
angle around — e z . Initial precession of Si induces pre- 
cession of S2 - mainly via the RKKY coupling. Gener- 
ally, response to the exchange field is slower than current- 
induced dynamics. Therefore, a difference in precession 
phase of S2 and Si appears, and configuration of SyAF 
deviates from the initial antiparallel one. This in turn 
enhances the STT acting on F2, which tends to switch 
S2. Its amplitude, however, is small in comparison to the 
strong RKKY coupling. Further scenario of the dynam- 
ics depends then on the external magnetic field. When 
77 app < 77 [Figs. UJb - e)] the Zeeman energy of S2 
has a maximum in the initial state and external mag- 
netic field tends to switch S2 to the opposite orienta- 
tion. Competition between the torques acting on SyAF 
results in out-of-plane precessions of both spins. After 
several precessions Si reaches the opposite static state, 
which is stable due to STT. In turn, S2 is only slightly 
affected by STT, and its dynamics is damped in the ex- 
ternal magnetic and RKKY exchange fields. In contrast, 
when 77 app > 77o [Figs. 0|f - i)] , Zeeman energy of F2 has 
a local minimum in the initial state, which stabilizes S2. 
Therefore, in a certain range of current density, SyAF 
does not switch but remains in self-sustained coherent in- 
plane precessions (red area in the upper part of Fig.[4ja)). 
For a sufficient current density, the SyAF structure be- 
comes destabilized and the precessional angle increases 
until the spins pass the (x, y)-plane. Consequently, the 
precessional angle decreases and spins of the SyAF are 
stabilized in the opposite state (t4-)- Moreover, as shown 
in Fig.HJc), the switching process for 77 app < 77 is con- 
nected with high distortion of SyAF configuration, where 
m in a certain point reaches its maximum value (corre- 
sponding to parallel orientation of both spins). Contrary, 
the m remains small for 77 app > 77 [Figs. HKg)] , and the 
effective magnetic moment of the free layer stays smaller 
than magnetic moment of a single layer. This might play 
an important role in applications of spin-torque devices 
based on CFLs. 

The two switching mechanisms described above dom- 
inate the current-induced dynamics when the current 
density is close to the dynamics threshold. For higher 
current densities, the nonlincarities in SyAF dynamics 
become more pronounced, which results in bistable be- 
havior of the dynamics, especially for 77 app < 77o and 
7 > 10 8 Acm~ 2 . In that region, the number of out-of- 
plane precessions before SyAF switching increases with 
the current density. However, their precessional angle 
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increases in time and consequently S\ might reach an 
out-of-plane static point slightly tilted away from the e x 
direction while S2 = e z remains in the layers plane. The 
out-of-plane static states (marked as ■(— f) have small re- 
sistance and appear as dark red spots in the diagram 
shown in Fig. Eta). 

In addition, from the analysis of the dispersion of pillar 
resistance (not shown) one finds that except of a narrow 
region close to the dynamics threshold with persistent 
in-plane precessions, no significant steady-state dynam- 
ics of SyAF element appears. As will be shown below, 
such a behavior might be observed when CFL becomes 
asymmetric (SyF free layer). 



C. Synthetic ferrimagnet 

Let us study now spin valve with SyF as a free layer, as- 
suming d\ = 4 nm and (I2 = 2 nm. Accordingly, -Hrkky 2 
remains 2 kOe while Hrkky 1 is reduced to 1 kOc. As in 
the case of SyAF, from the averaged time-dependent part 
of the pillar resistance we have constructed a diagram 
presenting current-induced dynamics, see Fig. [SJa). The 
static part of resistance is now R sp = 19.80 ffhn 2 . The 
diagram has some features similar to those studied in 
the previous subsection. However, the maximum critical 
current is shifted towards negative values of H &vv , even 
if magnetostatic interaction between magnetic layers is 
neglected. This asymmetry is caused by the difference 
in exchange and demagnetization fields acting on layers 
Fi and F2. Moreover, this difference leads to more com- 
plex dynamics of the CFL's spins than that in the case 
of SyAF. 

Generally, there are several dynamic regimes to be dis- 
tinguished. The first one is the region of switching from 
J.j' configuration to the opposite one, XI, which is located 
at largest values of H app in the diagram. Mechanism of 
the switching is similar to that of SyAF shown in Figs.|3|f 
- i), where CFL changes its configuration just via in-plane 
precessional states with a small value of m (weak distor- 
tion of the antiparallel alignment of Si and S^)- Further- 
more, the darker area above Hq indicates one of the pos- 
sible self-sustained dynamic regimes of SyF, i.e. the in- 
plane precessions (IPP); see Figs. [5fb - e). This preces- 
sional regime starts directly after the SyF switching, and 
Si and S2 precess around e z and — e z , respectively. Due 
to different effective fields in Fi and F 2 , and energy gains 
due to STT, the spins precess with different precessional 
angles [Fig. \5fc)] and consequently different frequencies. 
Because of the strong interlayer coupling and spin trans- 
fer between the layers, amplitudes of their precessions 
are periodically modulated in time. This modulation ap- 
pears also in the time dependence of pillar resistance. 
Conversely, below H the dynamic is dominated by large 
angle out-of-plane precessions (OPP) of both spins, as 
shown in Figs.^f - i). This dynamic state is connected 
with a strong distortion of the antiparallel CFL config- 
uration, i.e. large value of m, and large variation of the 
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FIG. 5. (color online) (a) Averaged resistance of Cu- 
IrMn(10)/Py(8)/Cu(8)/Co(4)/Cu(l)/Co(2)-Cu spin valve 
pillar with a SyF free layer, presented as a function of cur- 
rent density and applied magnetic field. Examples of current- 
induced dynamics for / = —3 x 10 8 Am -1 and H app = 200 Oe 
(b - e) and H app = -400 Oe (f - i). Panels (b) and (f) show 
dynamics of z-components of Si and S2, (c) and (g) present 
the overal magnetization of the free layer, (d) and (h) show 
the corresponding variation of pillar resistance, (e) and (i) 
show spin trajectories of Si (red solid line) and S2 (black 
dashed line) taken from a time interval as large as 30 ns after 
100 ns of initial equilibration. 

resistance. From Fig. [5ji) one can see that trajectories 
of Si and S2 are rather complicated including both IPP 
and OPP regimes with dominant OPP. 

D. Power spectral density 

From the analysis of current-induced dynamics we 
found that self-sustained dynamics in structures with 
SyF free layer is much richer than that in systems with 
SyAF free layer [see Figs. [SJb - i)]. Therefore, in this 
section we restrict ourselves to dynamic regimes of the 
SyF free layer only. More specifically, we shall examine 
the power spectral density (PSD) as a function of current 
density and external magnetic field. 
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FIG. 6. Power spectral density calculated for the spin valve 
with SyF free layer at T e g = 5 K and H app = —400 Oe (a) and 
200 Oe (b). (c) and (d) show steady time evolution of spins z- 
components in a time window of 10 ns after the equilibration 
at H app = -400 Oe and I = -2.8 x 10 8 AcnT 2 . (e) and (f) 
show steady time evolution of the time-dependent part of spin 
valve resistance in a time window of 20 ns after the equilibra- 
tion at // app = -400 Oe and I = -3.6 x 10 s Acm~ 2 (both Si 
and S2 precess out of layer plane) and I = —3.8 x 10 8 Acm -2 
(S>2 performs out-of-plane precessions while Si precesses in 
the layer's plane), respectively. Panels (g) and (h) depict 
trajectories of Si (red solid line) and S2 (black dashed line) 
corresponding to resistance oscillations (e) and (f), respectiv- 
elly. 



In the simulation we started from 1 = and changed 
current density in steps AI = 10 6 Acm~ 2 at a fixed ap- 
plied field. As before, to protect the SyF dynamics from 
collapsing into collinear static state, we assumed small 
thermal fluctuations corresponding to T ft = 5K. At 
each step we simulated the dynamics of coupled CFL's 
spins and calculated PSD. As in Ref. H^, we assumed 
that the input current is split between a load with resis- 
tance i?L and pillar with resistance R sp + R(t). Hence, 
voltage on the pillar has been calculated as U(t) = 
IR{t)/[l + Rs P /(R L S)], where we assumed R L = 50fl, 
and S is the cross-section of the pillar (ellipsoid with the 
major and minor axes equal to 130 and 60 ran, respec- 
tively). Then, at a given / we calculated voltage in the 
frequency domain, U(f), using fast Fourier transforma- 
tion over the period £fft = 50 ns following the equili- 
bration time of t cq = 30 ns. The power spectral density 
has been defined as PSD(/) = 2U 2 {f)/{R L Af), where 

A/ = 1/tpFT. 

Figures HJa) and (b) show PSD calculated at iJ ap p = 



—400 Oe and 200 Oe, respectively. The former case cor- 
responds to that part of the diagram in Fig. [SJa), which 
includes OPP modes, while in the latter case we observed 
IPP only. Let us analyze first the situation in Fig. . 
When current passes through the corresponding thresh- 
old value, both spins start precessing in the layers' plane, 
similarly as shown in Fig. EJb). Apart from the main 
peak in the PSD at / ~ 40 GHz, two additional minor 
peaks close to / ~ 20 GHz are visible. We attribute 
them to the oscillations of precessional amplitudes of 
both spins. With increasing amplitude of the current 
density, the precessional angles of both spins increase and 
their precessional frequencies slightly decrease. More- 
over, with increasing current the frequencies of the minor 
peaks become closer and closer, until they finally coin- 
cide. At this point the PSD becomes widely distributed 
along the whole range of observed frequencies, which is an 
evidence of noisy variation of the resistance. An example 
of spin dynamics in this region is shown in Figs. |5{c) and 
(d) which have been taken in a time window as large as 



10 ns after the equilibration period for H a 



-400 Oe 



and / = —2.8 x 10 8 Acm [within the broad feature 
of PSD in Fig. |B{a)]. Firstly, the figures show that S2 
starts to perform out-of-plane precessions as a result of 
the competition between STT and RKKY coupling. Sec- 
ondly, one can note thermally activated random transi- 
tions of £?i between OPP and IPP modes. These ran- 
dom transitions modify OPP precessions of S2 as well. 
Simultaneous dynamics of both spins causes chaotic vari- 
ation of spin valve resistance and broadens the PSD. The 
quasi-chaotic feature of the spin dynamics in this range of 
current densities can be seen also on the spin trajectories, 
which cover almost the whole sphere (not shown). 

Further increase in current density leads to stabiliza- 
tion of the OPP mode of Si . Hence spin valve resistance 
becomes more periodic [see Fig. EJe)] and PSD reveals 
a narrow peak again. Since both spins perform rather 
complicated dynamics including IPP but dominated by 
OPP regime [see Fig. IHg)], we observe a blue-shift in 
PSD with current, which is connected with a decrease in 
the precessional angles. However, at a certain value of 
/ we notice an abrupt drop in the peak's frequency. At 
this current density the STT acting on the left interface 
of layer Fi starts to dominate the dynamics of Si and en- 
ables only small angle IPPs along the So direction, which 
modifies the trajectory of S>2- S2 still remains in the OPP 
regime [see Fig. |HJh)] and hence the blue-shift with cur- 
rent appears. The fact that IPP of Si still influence the 
dynamics of the whole SyF is also shown in Fig. [6jf), 
which presents the dynamic part of the spin valve resis- 
tance at I = -3.8 x 10 8 Acm~ 2 and H app = -400 Oe. As 
a result of IPPs of Si , amplitude of the resistance varies 
periodically. In addition, comparison of Figs. |H|c) and 
(f ) shows that the simultaneous OPPs of both spins lead 
to stronger variation of the resistance than in the case 
when the layers are in the IPP state. 



Contrary, at H a 



200 Oe one observes only IPP 



modes of both spins similar to those shown in Fig. [SJe) 
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The in-plane precessional angle increases with current 
density and hence the peak frequency in PSD decreases 
and becomes broader. In real systems, however, one 
might expect the peaks narrower than those obtained in 
the macrospin simulations, as observed in standard spin 
valves with a simple free laye r 30 ! 31 . 



IV. CRITICAL CURRENTS 

First, we derive some approximate expressions for crit- 
ical current density needed to induce dynamics of CFL, 
derived from linearized LLG equation. In metallic struc- 
tures, the out-of-plane torque components are generally 
much smaller than the in-plane ones, and therefore will 
be omitted in the analytical considerations of this section 

(6<°\6«tf> -►()). 

The coupled LLG equations in spherical coordinates 
can be then written as 



dS 
~dt 



1 



1 



M 



(6) 



where S = [6\, <f>i, 02, 4>2) T is a 4-dimensional column vec- 
tor which describes spin orientation in both layers consti- 
tuting the CFL, and v = (vie, i>i0, V28, V2<j>) T stands for 
the torque components, v t g = T t ■ e.^ and v,^ = T,i ■ e^, 

with e it p = (e z x Si)/ sin6>i and e ie = {Si x e^)/ sin6»j 
denoting unit vectors in local spherical coordinates asso- 
ciated with Si. In turn, the 4x4 matrix M takes the 
form 



M 




(7) 



Static points of the CFL dynamics have to satisfy Vig = 
and Vi$ - for both i = 1 and i = 2. These conditions 
are obeyed in all collinear configurations, i.e. = 0,7r. 
Additional four trivial static points can be found in the 
out-of-plane configurations with 8i = tt/2 and 4>i = 0, 7r. 

Following Ref. [2J we linearize Eq. (J6j> by expanding v 
into a series around the static points, which leads to 



dS 



1 



1 



M J 6v 



(8) 



where J is a Jacobian matrix of dvi/dSj components. 
The matrix product MJ defines here the dynamic matrix 
D = M ■ J. This matrix allows one to study stability of 
CFL's spins in their static points. If Tr {D} is negative, 
the static point is stable, otherwise it is unstable. Hence, 
the condition for critical current ifM- Tr {D} =0. 

To obtain threshold current for dynamics onset of in- 
dividual spins in the CFL, we assume first that one of 
the spins is fixed in its initial position and investigate 
stability of the second spin. The dynamic matrix D be- 
comes then reduced to a 2 x 2 matrix. Considering initial 



position of SAF with Si — — S2 = — e z (i.e. 9\ — n and 

62 = 0), marked as 4<t, and polarizer So = e 2 , the stabil- 

z-4t 



ity condition leads to the following critical currents J^i 

■at 

c 2 

H M s di 



and I„ for S\ and S2, respectively: 



-it 



-,(2) 



— -^ext + -^ani — ^RKKY 1 + H° 



with H^ t = H^p - Hf - H 21 , and 



(9) 



r-lt 



fi M s d 2 



H, 



2# 



f H anl — -ffRKKY 2 + H 2 

(10) 

Hi 2 . In both above expres- 

and a^ 1 ' are taken in the considered static 
point, while the demagnetization field for the i-th layer 
is given by 



with H^l = H app 

(0) (2) 

sions a\ , a\ 



±J - Z 



Hf = 



H 



in 



— Hi, 



(11) 



Analogically, one can derive similar formulas for critical 
currents in the opposite (tj.) magnetic configuration of 
the CFL. 

Now we relax the assumption that one of the spins is 
fixed, and consider both spins of the CFL as free. Then, 
we calculate the trace of the whole 4x4 matrix, which 
leads to the following expression for critical current desta- 
bilizing the whole CFL: 



r 4-t 

'cCFL 



fioM s did 2 



d 2 (<4 0) +a¥>) + d 1 a. 



H 



2H„ 



— ^RKKY 1 — ffRKKY 2 + #1 + H 2 



[12) 



where H^. = Hf - Hf + H 21 + Hi 2 . Since the spins 
of CFL are antiparallel in the considered static point, 
-^cCfl 1S independent of external magnetic field. The 
above equation describes the critical current at which 
the CFL is destabilized as a rigid structure (unaffected 
by external magnetic field along the z-axis). 

Numerical calculations presented below show that crit- 
ical current is usually smaller than that given by Eq.(12). 
Apparently, as shown by numerical simulations, there is 
a phase shift in initial precessions of S\ and S2. Such a 
phase shift slightly perturbs initial antiparallel configu- 
ration and might reduce the critical current for the dy- 
namics onset. 

Similar formula also holds for the opposite configura- 
tion (ti), where the critical current is given by 



T 1i 

z cCFL 



fj, M s did2 



-H, 



n 



d 2 {af ) -af ) )-d 1 a ( i ) 

^RKKY 1 — ffRKKY 2 + #1 + H% 



2iT a 



(13) 
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with = Hf - Hf - Ef Hf. 

Now we discuss the theoretical results on critical cur- 
rents in the context of those following from numerical 
simulations. Let us consider first the critical currents for 
individual spins of the SyAF free layer, assuming that the 
second spin remains stable in its initial position, Eqs. (|9]) 
and (|10p . The corresponding results obtained from the 
formula derived above are presented in Table U where 
we have omitted a weak dependence on H &pp . For the 
studied structure with SyAF free layer, if is negative 
while 1*2 is positive. From our analysis follows, that 
the current density at which dynamics appears in the 
simulations (Fig. Eta)) is higher than that given by if. 
However, we checked numerically that the critical value 
I*i agrees with the critical current obtained from sim- 
ulations when assuming Si as free and fixing S-2 along 





SyAF 

It 


tl 


SyF 

It 


ti 


Id 


-0.31 


0.43 


-0.54 


0.63 


Ic2 


0.98 


-0.46 


2.33 


-0.78 


IcCFL 


-0.86 


0.87 


-2.37 


-9.18 



TABLE I. Critical current densities in the units of 10 Acm - , 
calculated according to equations (|10p . and (|12p for both 
SyAF and SyF free layers. 

Following the above discussion of the CFL dynamics, 
one can understand the shift of the threshold current 
as follows. Initially, when the current density exceeds 
if , Si becomes destabilized. Then, S2 responses to the 
initial dynamics of Si with similar coherent precession. 
However, S2 should still be stable in its initial position 
at this current density and common precessions of the 
two coupled spin moments damps the initial dynamics. 
Accordingly, SyAF ends up in the closest static state 
(It)- However, as the current density increases, the ini- 
tial precessions of Si become more pronounced, which 
in turn means that the initial antiparallel configuration 
becomes distorted and S2 becomes destabilized. This re- 
sults in coupled dynamics of both spins and finally leads 
to switching of the SyAF structure. 

On the other hand, we have also calculated the crit- 
ical current for the whole SyAF structure according to 
Eq. (TT2|) . and for the given structure we got lf CFL shown 
in Fig. Sta) by the dashed vertical line (see also Table 
I). Equation (jT!?j) describes stability of the whole CFL, 
and since the interlayer coupling is strong, lf CFL corre- 
sponds to the current density at which both spins become 
destabilized simultaneously preserving their antiparallel 
orientation. As can be seen in Fig. HJa) , this is the max- 
imum threshold current density for current-induced dy- 
namics. Because the rigid structure consisting of two 
antiparallel spins is not influenced by an external homo- 
geneous magnetic field, there is no dependence of ifcFL 



on -ffapp- Nevertheless, from our numerical simulation 
follows that the threshold current for the SyAF dynam- 
ics, I tbx , obeys the condition \lf\ < \I thl \ < \lf CFL \, 

provided that \lf\ < \lf\ or ( as m our case) if has 
different sign. 

When the SyAF is in the tl configuration, the spin 
accumulation and spin current are different from those 
in the |t configuration (at the same voltage). This in 
turn leads to different spin torques, which is the rea- 
son why the critical currents destabilizing tl state are 
different from those for ||~, as shown in Table |TJ From 
the critical currents one can expect relatively symmetric 
hysteresis with applied current in structures with SyAF. 
In contrast, I^ FL for the SyF is negative, similarly as 
lfc FL , k 11 ^ ^ is significantly larger, which indicates lack 
of hysteresis. To compare switching of the SyAF and 
SyF free layers from the |t to tl configurations with the 
opposite one (tl to |t), we have simulated dynamics of 
the corresponding CFLs assuming H &pp = and varying 
current density. The simulations have been performed 
in the quasistatic regime, i.e., for each value of current 
density the spin dynamics was first equilibrated for 50 ns 
and then averaged values of spin components and pillar 
resistance were calculated from the data taken for the 
next 30 ns of dynamics. In order to prevent the system 
from collapsing into a static state with zero torque, we 
have included a thermal stochastic field corresponding to 
T ff = 5K [see Eq. (|3J)]. Starting from / = and go- 
ing first towards negative currents we have constructed 
the current dependence of the averaged resistance and 
related z-components of both spins, as shown in Fig. [7] 
For both SyAF [Figs. [TJa -c)] and SyF [Figs, [^d - f)] 
free layers, one can see relatively symmetric hysteresis 
with the current density. In both cases direct switching 
from |t to tl state occurs at a current density comparable 
to if cfl- I n contrast, in the case of SyF free layer, the 
second transition (|| to |t) appears at a current density 
which is very different from that predicted by the lin- 
earized LLG model. Moreover, in both cases switching 
from tl to |t state does not appear directly, but through 
some precessional states. More precisely, as the positive 
current density increases, both spins start precessing in 
the layers' plane prior to switching. The in-plane preces- 
sions are connected with a significant drop in the resis- 
tance and with a reduction of the s 2 -components. The 
range of IPP regime is particularly large in the case of 
SyF. From the analysis of spins' trajectories one may 
conclude that the angle of IPPs increases with increasing 
current density, and after exceeding a certain threshold 
angle CFL switches to the |t configuration. 

The other factor giving rise to the the difference in 
switching from tl to |t and from |t to tl follows from the 
fact that the magnetostatic interaction of the CFL's lay- 
ers with the polarizer is different in the |t and tl states. 
To prove this we have constructed analogical hysteresis 
loops for SyAF and SyF free layers disregarding magne- 
tostatic interaction with the Fq layer; see Figs. [TJg) and 
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FIG. 7. Hysteresis loops of the resistance for the studied 
pillars with SyAF (a) and SyF (d) free layers. Panel (b) and 
(c) depict spin dynamics of Si and S2 in SyAF, respectivelly, 
corresponding to resistance loop (a). Panels (e) and (f) show 
dynamics of Si and S2 in SyF, respectivelly, corresponding 
to resistance loop (d). The initial point of each hysteresis 
loop is marked with a dot. The arrows indicate direction of 
the current change. Figures (g) and (h) correspond to the 
upper parts of (a) and (d), in which however the effects due 
to magnetostatic field of the reference layer to the CFL spins 
have been omitted. 



(h) . For both SyAF and SyF free layers we observe now 
large decrease in R for both switchings. This implies that 
both switchings are realized via in-plane precessions, in 
contrast to the case when Fo influences the CFL dynam- 
ics via the corresponding magnetostatic field. While the 
hysteresis loop for SyAF remains symmetric, the one for 
SyF becomes highly asymmetric. The asymmetry of SyF 
loop is due to a significant asymmetry of STT in and 
4-t states, which was previously shaded by the magneto- 
static coupling with the layer Fq. 



V. DISCUSSION AND CONCLUSIONS 

We have studied current-induced dynamics of SyAF 
and SyF composite free layers. By means of numerical 
simulations we identified variety of dynamical regimes. 
The most significant difference between dynamics of 
SyAF and SyF free layers concerns the evidence of self- 
sustained dynamics of both CFL spins. While in the 
case of SyAF only coupled in-plane precessions in a nar- 
row window of external parameters (H app and I) are ob- 
served, SyF free layer reveals more complex and richer 
dynamics, with the possibility of coupled out-of-plane 
precessions which might be interesting from the applica- 



tion point of view. Furthermore, as shown by numerical 
simulations, both SyAF and SyF are switchable back and 
forth without the need of external magnetic field. For 
SyAF element two possible ways of switching have been 
identified. Since they lead to different switching times, 
their identification might be crucial for optimization of 
switching in real devices with SyAF free layers. However, 
one has to note that the diagrams shown in Figs. 2] and [5] 
may be changed when magnetization in CFL becomes 
non-homogeneous. 

A disadvantage of the studied structures is their rel- 
atively low efficiency of switching, i.e. high amplitude 
of critical current and long switching time. In order to 
show more sophisticated ways of tuning the CFL devices, 
we have analyzed critical currents derived from the lin- 
earized LLC equation. The formula (|12[) has been iden- 
tified as the maximum value of critical current at which 
dynamics of the CFL structure should be observed. This 
formula reveals some basic dependence of critical current 
on spin valve parameters, and therefore might be use- 
ful as an initial tool for its tuning. However, in some 
cases non-linear effects in CFL dynamics might com- 
pletely change the process of CFL switching, as shown 
by the presented numerical simulations. But the effects 
of non-linear dynamics go beyond the simple approach of 
linearized LLC equation, and their study requires more 
sophisticated nontrivial methods and/or numerical sim- 
ulations. 
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Appendix A: Transformations of spin current 

Torque acting on the left interface of Fi is calculated 
from x and y components of j' x = T(0i, fa).ji, where ji is 
spin current vector in Ni layer (written in global frame; 
shown in Fig. Q}, and T(0i, fa) = R^ (-0i)R* {fa - 7r/2), 
where R 9 (a) is the matrix of rotation by angle a along 
the axis q in the counterclockwise direction when look- 
ing towards origin of the coordinate system. Hence 
components can be written as 

fix = hx sin fa - ji y cos fa , (Ala) 
fly = (jix cos fa + ji y sin fa ) cos X - j lx sin 0i , (Alb) 
fiz = ihx cos fa + jiy sin fa ) sin 0i + j lz cos 0i , (Ale) 

where (9i,fa) are spherical coordinates of Si in the 
global frame. Similarly, we define torques' amplitudes 
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on the left interface of F2 from the components of trans- 
formed spin current vector j 2 " = T(# 2 , 02) 02- In this 
case, however, J2 is not written in the global frame, but 
in the local coordination system coordinate system con- 
nected with Si. To rotate local coordinate system of 
Si to local coordinate system of S2 we need to know 
spherical angles 62 and 02 of vector S2 in the local coor- 
dinate system of Si. This might be done by transform- 
ing first S2 vector to local coordinate system of Si as 
S' 2 = T(#i,0i) • S 2 an d calculate its angles 9 2 and 2 . 
Then we can calculate components of j 2 " similarly as for 



the left interface 

j'i'x = hx sin <f> 2 - hy cos (j> 2 , (A2a) 
j'2'y = {hx cos (j) 2 + hy sin 02 ) cos 6» 2 - j 2z sin 2 , ( A2b) 
j'2'z = (hx cos 4>2 + hy sin 02 ) sin 9 2 + hz cos 9 2 , ( A2c) 

Equation in N 2 , which is adjacent non-magnetic inter- 
face from the right-hand side of Fi, are written in local 
coordinate system of Si. To apply the definition of ai2 
and 612 we need to rotate the local coordinate system so, 
that its y-axis will lie in the layer given by vectors Si 
and S'2. This might be done by single rotation of local 
coordinate system around its z-axis by angle 02 — t/2, 
j 2 ' = R z (02 - 7r/2).j 2 , where 

j'L = hx sin 02 - hy cos 02 , (A3a) 
j'iy = hx cos 02 + :)2 V sin 2 , (A3b) 
3'L=hz- (A3c) 

Note, angle 02 is calculated for vector S2 transformed 
into coordinate system of Si as in previous case. 
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